hep-ph/931020^ 



BI-TP 93/52 



GAUGE-INDEPENDENT EXTRACTION OF 
THE NEXT-TO-LEADING-ORDER DEBYE MASS 
FROM THE GLUON PROPAGATOR 



A. K. RebhanQ 

Fakultat fur Physik, Universitat Bielefeld 
D-33501 Bielefeld, Germany 



ABSTRACT 

It is shown that by defining the Debye mass through the relevant pole 
of the static gluon propagator rather than the zero-momentum limit of 
the time-time component of the gluon self-energy, a gauge- independent 
result for the next-to-leading order correction can be derived upon resum- 
mation of hard-thermal-loop contributions. The result turns out to be 
logarithmically sensitive to the magnetic screening mass. 
[Talk given at the 3rd Workshop on Thermal Field Theories and their 
Applications, August 15-27, 1993, Banff, Alberta, Canada] 



October 1993 



1 On leave of absence from Institut fiir Theoretische Physik der Technischen Universitat Wien, A- 
1040 Vienna, Austria. Address after November 1, 1993: DESY, Gruppe Theorie, D-22603 Hamburg, 
Germany. 



GAUGE-INDEPENDENT EXTRACTION OF 
THE NEXT-TO-LEADING-ORDER DEBYE MASS 
FROM THE GLUON PROPAGATOR 



A. K. REBHAN 

Fakultat fur Physik, Universitdt Bielefeld 
D-33501 Bielefeld, Germany 



ABSTRACT 

It is shown that by defining the Debye mass through the relevant pole of the static 
gluon propagator rather than the zero-momentum limit of the time-time component 
of the gluon self-energy, a gauge-independent result for the next-to-leading order 
correction can be derived upon resummation of hard-thermal- loop contributions. 
The result turns out to be logarithmically sensitive to the magnetic screening mass. 



1. The Hot QCD Debye Mass Puzzle 

In QED, the electric permittivity e(u, k) is given by the time-time component 
of the photon self-energy, 

/ t \ i . n 00(^ fc ) / -I \ 

e{u,k) = l + — 2 • (1) 

In the static limit to — 0, this is the screening factor of longitudinal electric fieldsB, 

<^W> = - Wr/m V\ (2) 

k 2 + n o(o,fc) 

where J M is a weak, conserved external current. The long-wavelength limit of 
Hoo(0, k) is usually identified with the Debye screening mass, m 2 el = f n o(0, k — > 0). 
In fact, in the high-temperature limit (T ^> k) the leading contribution 

e 2 T 2 

n 00 (o,fc) = ^- (3) 

provides a momentum-independent mass term, which modifies the classical Coulomb 
potential by a factor e _mcl r . The leading high-temperature result in QCD is very 
similar; one just has to replace e 2 — > g 2 (N + Nf/2) for color group SU(A^) and Nf 
fermionsEJ. 

In QED, higher-order corrections to the long-wavelength limit of n o(0, k) are 
known through an exact relation! to the equation of state 

n 00 (0,A:^0)=e 2 g2P ^ r) , (4) 



where P(fi, T) is the pressure and \i the chemical potential. Since the first three 
terms in the perturbation expansion of P are known, one also knows the first three 



terms in (g). With fj, = OjB 



e 2 T 2 / 3e 2 Vie 
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n»(o,*-o) = — + + (5 ) 

In QCD, one would expect the first correction term to be of relative order g 
rather than g 2 , because of the "plasmon effect"H. However, after ring resummation 
a gauge-fixing dependent results is found for n o(0, k — > 0) in covariant gauges 



n 00 (0, k^ 0) = m 2 0cl 1 + a —J—— g (6) 




with gauge fixing parameter a, which shows that it cannot be identified with a 
directly measurable quantity that the Debye mass presumably should be. 

The reason for the gauge dependence of (||) of course is that the self-energy of 
non-Abelian gauge fields is a gauge variant quantity. In QED, on the other hand, 
the self-energy of photons is directly related to the correlator of the gauge-invariant 
electromagnetic field strengths, so this problem does not arise there. 

It was therefore arguedi'lS that the non-Abelian Debye mass cannot in general 
be derived from the gluon propagator except in the temporal axial gauge, where the 
time-time component of the propagator is again directly related to the correlator of 
two chromoelectric field operators. 

The temporal axial gauge, however, is notoriously difficult at finite temperature^, 
and actually the principal-value prescription commonly used to deal with its singu- 
larities at zero frequencies has proved to be flawedQ. These singularities moreover 
prevent a straightforward implementation of ring resummation, because the static 
mode cannot be isolated. A first attempta of a ring-resummed calculation of the 
non-Abelian Debye mass yielded a negative correction term to the Debye mass of 
relative order g, but this could not be reproduced by taking the temporal limit of 
a corresponding calculation in general axial gauged, which by the way gave a posi- 
tive result. The former was consequently withdrawn and replaced! by a calculation 
which resums the asymptotic gluon mass rather than the leading-order Debye mass 
(not taking into account vertex corrections, however) . This yielded again a negative 
correction to the Debye mass as defined through the infrared limit of the self-energy 
in temporal axial gauge, to wit, 



lC G (0, k - 0) « ml, (l - J-^(iV + N f /2)j . (7) 

Leaving aside for the moment the open questions about pole prescriptions in 
temporal axial gauge and even which resummation scheme should be employed, 
there still remains the question whether the analysis of correlators of chromoelectric 
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field operators really guarantees to give gauge independent answers where the gluon 
propagator evidently failed. The temporal axial gauge is singled out as the one 
where no higher vertex functions are needed to obtain this correlator, but one could 
of course use any gauge. Because the chromoelectric field is not a gauge-invariant 
operator, there is a priori no reason to expect gauge fixing independence of its 
correlation functions. Indeed, one finds that under a change of gauge condition 
f^A^ — > {f/jjA- Sf^A^ the correlator of two chromoelectric field operators varies 
according tea 

8(E«{x)Et{y)) = -gf abc J d i z(E b J (x)c c (x)c d (z)6f ll A d ^z)E e k (y)} + {a,j,x~ e } k,y) } 

(8) 

where c and c are Faddeev-Popov ghost fields, which make their appearance even 
in gauges which are otherwise ghost-free. 

Thus one is not necessarily on safe grounds by studying the correlation of field 
operators rather than of the gauge fields. On the other hand, if one succeeded in 
extracting gauge independent information from the gauge variant quantity n o, then 
this information should equally be found in the correlator of electric field operators, 
since in a particular gauge (the temporal one) the two are directly related. 



2. Changing the Definition of the Debye Mass 

A strong hint that the very definition of the Debye mass as n o(0, k — > 0) might 
not be sufficient beyond leading order is seen by the consequences of keeping the 
term proportional to k 2 in the resummed result (B). Then£9 



Mo,* - o) = < L (i + a^^ s | - wi/svTlvy^ + " " (9) 

which, when inserted into (^j) leads to a different mass term besides an over-all 
factor that is constant. This does not remove the gauge dependence of the putative 
correction to the Debye mass (it just replaces a by a + |)0, but it does tell that 
the fc-dependence of I1 00 (0, k) still can change things! 

Let us therefore go back to the linear response formula for longitudinal electric 
fields, eq. (|^). Actually, this formula is validE-H also in the non-Abelian case, since 
with a single source J there is also only a single direction in color space to which 
the gauge potentials can point, and the nonlinear terms in the chromoelectric field 
strength vanish trivially. If Jo is a static point charge Q located at the origin, then 
the Fourier transform of eq. @ equals minus the gradient of a potential $ given by 



$(r) = Qj 



d 3 k 



e 



kr 



(27r) 3 k 2 + Tl 00 (k = 0, k) 

Q ,oo e ikr _ e -ikr fc & 



(2tt) 2 7-oo 2ir k 2 + U 00 (0,k)' 



(10) 
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The usual definition of the Debye mass as vr? d = f U 00 (0,k — ► 0) is commonly 
motivated by saying that when r is very large, the dominant contribution to the 
integral (|IUp comes from k = 0. This is not quite correct, however. Inserting for 
instance the leading high-temperature result (|3|), which is a simple constant mass 
squared, one can readily evaluate (|10|) by appropriately closing the contour in the 
complex fc-plane wherein the integrand has simple poles at k = ±zmo e i., where only 
Re/c = 0. If II 00 (0, k) depends on k, as it will be the case in general (i.e., beyond 
leading order) , then this dependence will be important to determine the location of 
the pole and thus the magnitude of the mass that will appear in the exponent of 
e~ mr . It is therefore rather obvious that one should define the Debye mass, which 
certainly should account for the actual exponential fall-off, self-consistently by the 
zeros of the denominator in (|TC)| ). I therefore propose the implicit definition 

m e \ = n o(o,fc)| fc2= _ m , (11) 

for the (chromo) electric Debye mass. 

The Debye mass is thus defined through the singularity of the propagator ap- 
pearing in (|T0| ) , and this propagator is in fact one of the components of the (static) 
gauge-field propagator for which formal arguments showing gauge independence of 
its poles can be derived from gauge fixing identities^. In non-Abelian gauge the- 
ories, the gauge dependence of the usual definition of the Debye mass through the 
zero-momentum limit of the self-energy is only to be expected because the leading 
high-temperature corrections move the pole away from strictly k = 0. 

One might now wonder why there appeared to be no problem with the old 
definition in the case of QED. In QED there is no problem of gauge dependence 
for the photon self-energy, for the latter is gauge independent through all orders of 
perturbation theory. Nevertheless, if one wants the Debye mass to truly describe 
the exponential fall-off of the electrostatic potential, then it is clear that one has to 
adopt the self-consistent definition (|TT|). But then there is a correction term to be 
added to (|5|) according to 



el. 



m e \ = n 00 (o, k - o) + [n 00 (o, k)\ k2= _ ml - n 00 (o, k - o) , (12) 

and this correction is 



m 2 el - n 00 (0, k -> 0) = m, 



2e 2 e 2 r 



Ocl. 



9ir 2 6tt 2 



In + lE 



+ 0(e 4 ) , (13) 



where jl is the mass scale introduced by dimensional regularization in which minimal 
subtraction has been performed. This correction term shows what has been missing 
in the first place when trying to identify ([5]) with a physical quantity. As a physical 
quantity, it has to be a renormalization-group invariant, which is indeed what ( |i~3|) 
brings about: the coefficient of the logarithmic term in (|i~3|) is exactly such that 
dml/dfl = because de/d{hxjl) = (3{e) = e 3 /(12n 2 ) + 0(e 5 ). 
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In fact, the importance of the fc-dependence of IIoo(0, k) for Debye screening 
has been recognized previously in the case of a degenerate electron gas, giving rise 
to the so-called Friedel oscillations". Yet the self-consistent definition ( |TT]) of the 
Debye mass has not been adopted before as far as I know. 



3. Next-to-Leading Order Calculation of the Non-Abelian Debye Mass 

A complete calculation of perturbative corrections to the high-temperature dis- 
persion laws beyond those determined by the gauge-independent hard thermal loops 
(HTL's) has been shown by Braaten and Pisarskili^ to require the resummation 
of all of the HTL contributions to self-energy and vertices. This rather involved 
resummation schemehas been applied in particular to determine the damping of 
collective excitations!^, and most recently also to compute the next-to-leading order 
contribution to the chromo dynamical plasma frequency t3. 

The problem of calculating the next-to-leading order contribution to the non- 
Abelian Debye mass is in fact a problem of the same kind. With the definition ( fjT|) , 
the Debye mass can be regarded as being given by the position of the plasmon pole 
when the frequency is lowered below the plasma frequency and eventually put to 
zero, whilst the wavevector becomes imaginary. 

Fortunately, with zero external frequencies the Braaten-Pisarski scheme can be 
simplified tremendouslycZl. In the imaginary time formalism, one may separate the 
modes of a resummed bosonic propagator l/[(27rnT) 2 + k 2 + II] into the static mode 
n = and the non-static ones, where the latter are automatically hard in the termi- 
nology of Braaten and Pisarski. Hence, only the static mode needs resummation of 
the HTL, which in the propagator is just the leading-order Debye mass m 0e i.. This 
separation of course makes analytic continuation to non-zero external frequencies 
rather impossible, but this is no shortcoming in the static case. With all the exter- 
nal frequencies being zero, all potentially soft lines of a diagram are static, and so 
are the vertices that would need resummation of HTL contributions. However, the 
latter vanish entirely in the purely static limit. 

Hence, as far as the relative order g correction to the static gluon self-energy 
is concerned, Braaten-Pisarski resummation boils down to a conventionally ring- 
resummedtl one-loop calculation. The static ring-resummed propagator in general 
covariant gauge reads 



p =0 



(14) 

Po=0 



and the complete next-to-leading order contribution to n o(0, k) is found as^l 



6 f d 3 ~ 2£ p f 1 1 



S n m(h = 0, k) = 9 mN^ WTWt J 



2 _|_ m 2 p2 



9 2 T 
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4m 2 -(k 2 + m 2 )[3 + 2pk/p 2 ] 2 , p 2 + 2pk ) 

+ jt^, T\ + a(k +m ) ^ k 15 

p (q + m^J p 4 (q z + rrr) J 

where m = m 0e i. and q = p + k. (Here dimensional regularization has been used 
when separating the static modes from the sum over Matsubara frequencies^; the 
limit e — > gives a regular expression because of the odd integration dimension.) 

The new definition (|TT|) requires to evaluate at k 2 = — m 2 . There the gauge 
parameter (a) dependent part vanishes algebraically, in accordance with the gauge 
fixing identities". However, closer inspection reveals that the integrals in ( |T5| ) 



develop "mass-shell" singularities, caused by massless transverse and massless un- 
physical modes in the gluon propagator. 

The third term of the integrand in (|I5|) is logarithmically singular as k 2 — > — m 2 , 
and the singularity is caused exclusively by the massless denominator in the spatially 



transverse part of the gluon propagator ([14]). A magnetic screening mass would 
remove this singularity, and because the latter is only logarithmic, the coefficient of 
the corresponding logarithm is determined by (|15|), 

att (n i\ g 2 Nm el T 2m el . 

<m 00 ((U) In- (16) 

up to terms that are regular as m magn — > 0. Assuming that m magn ~ gm^, the 
next-to-leading order contribution to m 2 x is found to be of order gln(l/g) rather 
than g, 

5ml N I 6 , 1 . 

and it is positive, at least at weak coupling g <C 1. Taken seriously for larger coupling 
g ro 1, the logarithm would eventually switch sign, but there the sublogarithmic 
terms would be of equal importance. 

Unfortunately, the sublogarithmic terms cannot be calculated completely, be- 
cause the presumed phenomenon of magnetic screening is nonperturbativeEj. How- 
ever, in order to obtain an estimate of those, let us assume that a simple replace- 
ment of 1/k 2 — > l/(k 2 + ) in the transverse part of the static propagator 



(|T4] ) correctly summarizes the effects at k ~ g 2 T. Then we may go on to evaluate 
the remaining contributions in (|I5|). Here one encounters a difficulty with the in- 
dependent term in Eq. (|15D, because by approaching the imaginary pole k 2 — > —m 2 , 
the explicit factor that apparently ensures gauge independence gets cancelled by a 
linear singularity in the momentum integral. Exactly the same phenomenon was 
encountered^ in the recalculation of plasmon damping rates in general covariant 
gauges. I have argued previously^ that this behaviour just reflects a singular, gauge 
dependent behaviour of the residue of the propagator rather than an actual gauge 
dependence of the pole determining the dispersion laws. Indeed, introducing an 
(unphysical) cut-off again moves the gauge dependence seemingly afflicting the pole 
position into the residue, while the correction to the pole position becomes indepen- 
dent of this infrared regularization. Alternatively, the gauge dependent contribu- 
tions can be avoided altogether by a quantization procedure that keeps unphysical 
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modes unthermalizedi. The gauge independent correction term then readsEl 

6 = ^ = 9 nJ J - L^- + \n2- 1 -)+ 0(g 2 ). (18) 
m 2 0el y y2N + N f 2ir\ m magn . 2 J {y ' K J 

In pure SU(2) lattice gaugetheory high-statistics results on Debye screening have 
rather recently been obtainedE3 at temperatures well above the critical temperature, 
finding a positive excess in the screening mass squared of 5 = +0.30(9). Inserting 
the parameters of this lattice calculation as well as an older lattice result for the 
magnetic screening massed, ( |T8|) yields 5 ~ +0.5, which comes remarkably close in 
view of g ~ 1. (The renormalized value of g used in this calculation appears to 
be in good agreement with independent lattice calculations of the SU(2) pressure 
when the latter is equated to the perturbative result .E3) 

In conclusion, by defining the Debye mass through the relevant pole of the 
static gluon propagator rather than the (gauge-dependent) zero-momentum limit 
of the gluon self-energy, a gauge-independent result at next-to-leading order has 
been derived after identifying and resumming the relevant hard-thermal-loop con- 
tributions. The location of the pole turns out to be logarithmically sensitive to 
the nonperturbative magnetic screening mass, but the coefficient of the correspond- 
ing logarithm can be calculated perturbatively. The latter is in fact related" to a 
logarithmic divergence ~ gm e \, ln(m e i r) encountered by Nadkarnit3 when trying to 
extract corrections to the non-Abelian Debye mass from the correlation function of 
two Polyakov loops. 

A similar logarithmic sensitivity to the scale g 2 T has been encountered in the 
calculations of plasmon and fermion damping^, which even appears in the Abelian 
case where there is no magnetic screening mass. In the case of the Debye mass, 
however, the origin of the logarithmic term is genuinely non-Abelian. Another dif- 
ference which is important in view of the discussions^ surrounding the calculations 
on damping is that the position of the pole defining the Debye mass is on the imag- 
inary axis of k, and it stays there when the corrections are included. This makes 
it rather unnatural to attempt anything else than a self-consistent procedure. As 
concerns the gauge dependences which in both calculations require regularization 
of mass-shell singularities (unless unphysical modes are frozeni), it should be kept 
in mind that they occur in the sublogarithmic terms which are strictly speaking 
beyond the reach of perturbation theory, at least in the non-Abelian case. It is 
therefore left open whether the ticklish "mass-shell singularities" caused by the 
massless unphysical modes in covariant gauges might not disappear in a calculation 
which is complete down to and including the order k ~ g 2 T. 

I would like to thank R. Baier, J. Kapusta, U. Kraemmer, R. Pisarski, K. 
Redlich, H. Schulz, and A. Smilga for useful discussions, and the organizers of 
the 3rd Workshop on Thermal Field Theories, in particular R. Kobes and G. Kun- 
statter, for this most enjoyable meeting. 
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